Neutron stars with isovector scalar correlations 



B. Liu^'^ H. M. Di Toro'^, V. Greco'''^ 

^ Center of Theoretical Nuclear Physics, National Laboratory of Heavy Ion Accelerator, Lanzhou 730000, China 
^ Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100039, China 
^ Department of Technical Physics, Peking University, Beijing 100871, China 
^ Laboratori Nazionali del Sud, Via S. Sofia 44, 
1-95123 Catania and University of Catania, Italy and 



5 



Cyclotron Institute, Texas A & M University, College Station, USA 



Neutron stars with the isovector scalar (S-field are studied in the framework of the relativistic mean 
field {RMF) approach in a pure nucleon plus lepton scheme. The 5-field leads to a larger repulsion 
in dense neutron-rich matter and to a definite splitting of proton and neutron effective masses. Both 
features are infiuencing the stability conditions of the neutron stars. Two parametrizations for the 
effective nonlinear Lagrangian density are used to calculate the nuclear equation of state (EOS) 
and the neutron star properties, and compared to correlated Dirac-Brueckner results. We conclude 
that in order to reproduce reasonable nuclear structure and neutron star properties within a RMF 
approach a density dependence of the coupling constants is required. 
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A Relativistic Mean Field (RMF) approach to nuclear matter with the coupling to an isovector 
scalar field, a virtual ao(980) (5- meson, has been studied for the asymmetric nuclear matter at low 
densities, including its linear response Q, and for heavy ion collisions at intermediate energies 

where larger density and momentum regions can be probed In this work we extend the analysis 

of the contribution of the ^-field in dense asymmetric matter to the influence on neutron star properties. 
We also want to test which effective interaction is more appropriate for the description of dense matter, 
including the symmetric part. In order to make a comparison, we use the two parametrizations for 
the effective nonlinear Lagrangian density to study the strongly isospin-asymmctric matter at high 
density regions. 

A Lagrangian density of the interacting many-particle system consisting of nucleons, isoscalar (scalar 
(T, vector cu), and isovector (scalar 6, vector p) mesons is the starting point of the RMF theory 



+ \{dj- d^8- m]5^) - \f,^F^'' - \g,.&\ (1) 

where (f) is the (/)-meson field, is the w-meson field, 6^ is p meson field, 5 is the isovector scalar 
field of the 5- meson. F^i, = d^uju — d^uj^, G^n = d^h^ ~ di,h^, and the U{4>) is a nonlinear potential 
of a meson : U{(j)) = ^acj)^ + jbcf)'*'. 

The field equations in a mean field approximation {MFA) are 



(«7^9^ - (M - g^c/) - gsTsSs) - g^l^uja - gpi°T^bo)ip 0, 
mlcj) + a(j)^ + bcjr' ^ g^ < ipip g^p^, 
mlujo = guj < ■07° >= du^P, 
"T-p^o = 5p < i'l^Tsip >= gpP3, 

"^^^3 = 55 < '0T-3V' >= gsPs3, (2) 

where pa = Pp — pn and p^a = Psp — Psn, P and ps are the baryon and the scalar densities, respectively. 
Neglecting the derivatives of mesons fields, the energy-momentum tensor in MFA is given by 

T^. = ihi^dA + [^m^c/.^ + U{<i>) + im^J^ „ ^ml^^u:^ - \m%b^]gp,. (3) 

The equation of state {EOS) for nuclear matter at T=0 is straightforwardly obtained from the 
energy-momentum tensor. The energy density has the form 

i—n,p 

and the pressure is 



where E* = aZ/s^ + M*"^ , i=P,n. The nucleon effective masses are, respectively 



M/ = M-g,0- 55,53, (6) 

and 

Mr^* ^ M - g„(t) + gsS3. (7) 
The nucleon chemical potentials pi are given in terms of the vector meson mean fields 



Pi — y kp^ + M* + gtoi^o T 9pba (+ proton, - neutron), (8) 

where the Fermi momentum fc^. of the nucleon is related to its density, kp^ = {Sn^ pi)^^^ . 

Since we are interested in the effects of the Nuclear Equation of State we will consider only pure 
nuclconic (+lcpton) neutron star structures^i.e. without strangeness bearing baryons and even de- 
confined quarks, see the recent nice review and the ref.Q. In particular we will use two models 
for the neutron star composition, pure neutron and /3-stable matter. In the latter case we limit the 
constituents to be neutrons, protons and electrons. Then the composition is determined by the request 



of charge neutrality and /3-equilibrium. The {n,p,e~) matter is indeed the most important /3-stable 
nucleon + lepton matter at low temperature. 

The chemical potential equilibrium condition for the {n,p, e^) system can be written as 



- • (9) 



The charge neutrality condition is 



Pe^ Pp^ Xpp , (10) 

where Xp = Z/A — Pp/p is the proton fraction (asymmetry parameter a = 1 — 2Xp), and p is the 
total baryon density. The electron density Pe in the ultrarelativistic limit for non-interacting electrons 
can be denoted as a function of its chemical potential 



Pe = ^t^l , (11) 



where p.e = \jkp^ + ml. The Xp can be obtained by using Eqs.(8), (9), (10) and (11). The Xp is 
related to the nuclear symmetry energy Esym 



Z^'pXp ~ [4£;,j,„.(p)(l - 2Xp)f - 0. (12) 

In presence of a coupling to an isovector-scalar 5-meson field, the expression for the symmetry energy 
has a simple transparent form, see 0, 0| : 

1 iP' 1 M* 



where M* = M - g^(j) and Ep* ^ \l k^, + M*^. The Esym and the EOS for the /3-stable {npe~) 
matter at T=0 can be estimated by using the obtained values of Xp. Equilibrium properties of the 
neutron stars can be finally studied by solving Tolmann-Oppenheimer-Volkov (TOV) equations [^ll(ll| 
inserting the derived nuclear EOS as an input, 7]. 

The isovector coupling constants, p-field and p + S cases, are fixed from the symmetry energy at 
saturation and from Dirac-Brueckner estimations, see the detailed discussions in refs.|3,y|- 

In order to make a comparison, two parameter sets for the isoscalar part are used. The first. Set A, 
is more suitable at high densities where it appears closer to various Dirac-Brueckner predictions. In 
fact recently this interaction has been used with success to describe reaction observables in RMF- 
transport simulations of relativistic heavy ion collisions, where high densities and momenta are reached 
mk TW s«o„d, S.t B. is take, from the iVL3 pa„„.etri.atL Q, obtained b.fl«„g p.op«.ie. 
of symmetric nuclear matter at saturation density and of finite nuclei. In ref.|^ it has been shown 
that the good description of finite nuclei, even exotic, is kept also when the isovector-scalar channel 
is included, normally not present in the NLi Lagrangian. 



The coupling constants, fi = gf /mf, i = a,uj,p,6, and the two parameters of the a self-interacting 
terms : A = a/g^ and B = b/g'^ are reported in Table 1. The corresponding properties of nuclear 
matter are listed in Table 2. 



Table 1. Parameter sets. 



Parameter 


Set A 


Set B 




NLp 


NLp(5 


NLp 




fa 


10.32924 


10.32924 


15.61225 


15.61225 




5.42341 


5.42341 


10.40068 


10.40068 




0.94999 


3.1500 


1.09659 


3.08509 


fs {fm') 


0.00 


2.500 


0.00 


2.400 


A ifm-^) 


0.03302 


0.03302 


0.00999 


0.00999 


B 


-0.00483 


-0.00483 


-0.002669 


-0.002669 



Table 2. Saturation properties of nuclear matter. 



Parameter sets 


A 


B 


pa ifm-^) 


0.16 


0.148 


E/A (MeV) 


-16.0 


-16.299 


K {MeV) 


240.0 


271.7 


Esy„r (MeV) 


31.3 


33.7 


M*/M 


0.75 


0.60 



We first use the two parametrizations to calculate the scalar and the vector potentials, and the 
binding energy E/A for symmetric nuclear matter (a=0.0) as a function of baryon density. The 
results are presented in Fig.l. 

The relativistic Dirac-Brueckncr-Hartrce-Fock [DBHF) approach is a microscopic model describing 
many-body system with correlations, that has been extensiyely used to study the nuclear matter 
properties [l2> ll^ llJj lid ll^ • In order to make a comparison, in Fig.l we also report the different 
results within the DBHF approaches, the relatiyistic Dirac-Brueckner calculations by Brockmann 
and Machleidt denoted as DBBM, the Dirac-Brueckner T-niatrix calculations Q|, denoted as 
DBT, and the Dirac-Brueckner results by ter Haar and Malfliet 

[l4 denoted as DBHM. 
From Fig. 1(a) we see that the scalar and vector potentials for the isoscalar channels, given by 
the Set B (i.e. the NL3 au! couplings) in low density regions are consistent with the correlated 
Dirac-Brueckner results, while the results given by the Set A are in better agreement in high density 
regions. 

The dotted line in Fig. 1(b) denotes the EOS of symmetric nuclear matter given by the Set B, 
in full overlap with the solid line given by the NLS interaction. It shows a nice agreement with 
correlated relativistic predictions at low densities but is clearly too repulsive with increasing density. 
At variance Fig. 1(b) also shows that, as expected, the EOS of symmetric nuclear matter given by 
the Set A is more consistent with that given by Dirac-Brueckner calculations, in particular for the 




FIG. 1: (a) Scalar and vector potentials vs the baryon density; (b) binding energy as a function of the baryon 
density for symmetric nuclear matter. See text. 



DBT and DBBM estimations, in high density regions up to about three-four times normal density. 
Such very different behaviors at high baryon density, in connection to the large difference in nuclcon 
effective masses, will strongly influence the neutron star structure. 

Relativistic heavy-ion collisions can provide crucial information about the EOS of nuclear matter. 
The investigation of nuclear EOS at high densities is one of driving force for study of heavy-ion 
reactions. The authors [3| use different DBHF approaches to study the collective flow of heavy-ion 
collisions. It is shown that the softer DBT choice is in better agreement with experimental data of 
relativistic collisions, at least up to a few AGeV beam energies, where densities up to 2.5po can be 
reached in the interacting zone. In general very accurate analyses of relativistic collisions data favor 
the predictions of a softer EOS at high densities [islllfll . We remark from Fig. 1(b) that for symmetric 
matter the DBT is quite close to our Set A parametrization, at least up to about S/Oq. 

For the isovector channels one can see from Eqs.(6) and (7) that the presence of the i5-field leads to 




FIG. 2: Neutron and proton effective masses vs. the baryon density for some values of the proton fraction: 
(a) by Set A and (b) by Set B, respectively. See text. 



proton and neutron effective mass spiitting. In Fig. 2 we present the the baryon density dependence of 
n, p effective masses for different proton fractions by the two parameter sets. The solid lines in Fig. 2 
are the nucleon effective mass for symmetric nuclear matter {Xp = 0.5). 

Fig. 2(a) shows that the proton and the neutron effective masses given by the Set A decrease slowly 
with increasing baryon density, at variance with the Set B case that presents a much faster decrease, 
Fig. 2(b). This main difference between the two, A and B, parametrizations, is actually coming from 
the isoscalar part. When coupled to the splitting due to the isovector (5-field it will have large effects 
on the n-star equilibrium features. 

The density dependence of symmetry energy for the two parameter sets is reported in Fig. 3. For 
both cases we see a similar behavior of Esym at sub-saturation densities for NLp and NLp6 models. 
With increasing baryon density p, however, the differences arising from the presence of the (5- meson 
in the isovector channel become more pronounced for both A and B, parametrizations. This is due to 




the quenching factor {M* /Ep)'^ for the attractive 5 contribution in Eq. (13) at high density regions, 
see refs.0,0]. 

We note that in spite of the same isovector couphng constants, at high density the Set B gives a 
larger symmetry energy. This is related to the larger contribution of the kinetic term in the r.h.s. 
of Eq.(13) due to the faster decrease of the effective nucleon mass. This represents a nice example 
of how the isovector part of the nuclear EOS can be influenced by the isoscalar channels due to the 
Fermi correlations. 

In the inserts of Fig. 3 we show the corresponding proton fractions Xp at /3-equilibrium, Eq. (12). 
Due to the stiffer nature of the symmetry energy in the NLpS cases in both parametrizations we see 
a decrease of the purca, i-e. of the baryon density corresponding to the value Xp = 1/9 that makes 
possible a direct Urea process, see 



At this point we can derive the predictions for equilibrium properties of neutron stars just solving 





FIG. 4: Mass of the neutron star as a function of the central density of the neutron star by Set A and Set B, 
respectively. 



the TOV equations |9|,ll0|. The main resuhs are presented in the Figs. 4 and 5. 

Fig. 4 displays the neutron-star mass as a function of the central density of the star given by the 
two parameter sets, for the two compositions, pure neutron and the /3-equilibrium (npe^) matter. 

Fig. 4(a) and (b) both show that the maximum masses of the /3-equilibrium star are smaller than in 
the pure neutron case due to the presence of a proton fraction and so of a smaller symmetry repulsion. 
Consistently the corresponding central densities are larger. A related effect is that the mass of the 
{npe~) star given by the NLpS model decreases more quickly with increasing density than that given 
by the NLp choice, i.e. we have a lower density instability onset. This is just because the introduction 
of the S coupling increases the equilibrium proton fraction at high baryon densities, see the inserts in 
Fig.3. 

The larger stiffness of the symmetry energy in the NLpS cases, in both parametrizations, can 
be directly seen in the fact that the corresponding curves are always above the ones without the S 



coupling. This implies larger maximum masses and smaller central densities. As a consequence we 
systematically see that the results for the (npe^) composition in the NLp6 models are approximately 
equivalent to the ones for the pure neutron matter in the NLp choices. 

Fig. 4(b) presents a qualitatively new feature of the Set B results: the lack of solution (maximum 
mass) in the /3-equilibrium case for the NLp5 model. The fast decrease of the neutron effective 
mass, see Fig. 2(b), prevents the chemical potential equilibrium condition for the {npe~) matter to be 
satisfied at densities around 3po- 
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FIG. 5: Mass of the neutron star as a function of the radius of the neutron star by Set A and Set B, respectively. 



Fig. 5 reports the correlation between neutron star mass and radius given by the two parameter sets, 
respectively, for the two cases, pure neutron and {npe~) matter. Fig. 5 shows that the contribution 
of the ^-field to the neutron stars at high density regions is quite remarkable. In particular we note 
that we systematically have larger masses and radii and lower central densities, as expected from the 
larger symmetry pressure. 

All the estimated maximum masses and the corresponding central densities and radii of the neutron 



Table 3. Maximum mass, corresponding radius and central density of the star given by the two 

parameter sets. 



Neutron star 


Properties 


Set A 


Set B 






NLp 


NLpS 


NLp 


NLp5 


pure neutron 


{Ms/MQ)maoo 


2.24 


2.56 


2.88 


3.12 




pc/pO 


5.75 


4.65 


4.12 


3.55 




R{km) 


11.44 


12.56 


13.70 


14.69 


(npe) matter 


{Ms/MQ)max 


2.14 


2.30 


2.83 






Pc/pO 


6.77 


6.49 


4.54 






R(km) 


10.80 


11.33 


13.25 





stars are reported in Table 3. 



We note that the difference between the two, A and parametrizations, in the neutron star pre- 
dictions, is largely due to the isoscalar structure of the interactions, see Fig.l. The B parametrization 
is much stiffer at high density regions and this leads to differences on the neutron star masses, radii 
and central densities. The comparison between the results given by the two sets shows that the set B 
(i.e. the NLZ forces) can be good at low densities below saturation density and it has serious prob- 
lems at high densities, while the set A is a good choice for EOS of nuclear matter at larger density 
regions, which is consistent with that pointed out by refs. l^^. In a sense all that just shows 
that we need density-dependent RMF parametrizations, the i3-type at low densities and A-type at 
high densities. This has been already emphasized in the work j2l|. While results of Density Depen- 
dent (DD) parametrizations and of the NL3 forces agree very well below the saturation density, the 
EOS of ZJZJ-interactions at supra-normal densities shows a much softer behavior, similar to DBHF 
calculations and in better agreement with Heavy Ion Collision [HIC) data, see also ref. |5|. 

In any case our calculations show that the 5-field provides significant contributions to the neutron 
star structure for the stiffness of the symmetry energy and the neutron/proton mass splitting in high 
density regions. The proton fraction in the /3-equilibrium matter is much larger than that in the no 
(5-field case. Moreover we can see from Fig. 4(b) and Fig. 5(b) that we do not have solutions in the 
(npe~) case for the set B with NLpd isovector interaction. This represents a quite dramatic effect of 
the splitting of neutron/proton effective masses at high densities, on top of the larger nucleon mass 
decrease of the B parametrizations, see Fig. 2. The neutron chemical potential is then not able to 
satisfy the /3-equilibrium conditions of the (npe~) matter. The contribution of the (5-field for strongly 
isospin-asymmetric dense matter is important and it cannot be neglected. 

With reference to neutron star properties we study the EOS for dense asymmetric matter in the 
RMF frame with two different parameter sets for the Lagrangian density. The set B is close to 
the NLi parametrization which has been proposed to describe finite nuclei properties. Since the 
proton and the neutron effective masses, especially the neutron effective mass, decrease quickly with 
increasing baryon density, the set B cannot provide the EOS needed in high density regions for the 
{npe~) star. This means that in general the B parametrization seems to have serious problems at high 
densities, as already remarked from HIC studies Though the B parametrization can be good 



for the EOS at low densities, below and around saturation, it is too stiff in high density regions, in 
particular there is no solution for the {npe~) case with the 5-field. So the set B is not suitable for the 
case of dense matter. We require a softer EOS at high densities, and indeed the softer Dirac-Brueckner 
predictions, in particular the DBT one, are in better agreement with relativistic collisions data. Our 
A parametrization, including the (5-isovector-scalar field, is quite close to the DBT and it has been 
shown to lead to good predictions in transport simulations for heavy ion collisions at intermediate 
energies P, IH 1^ - It appears then quite appropriate for the nucleonic part of the approach to neutron 
star properties. In this respect we note that quite extended neutron star structure calculations have 
been recently performed just using our Set A Lagrangian, [2^ . 
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